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SECTION  1 


INTRODUCTION 


Nuclear  exolosions  at  or  near  the  surface  of  the 
earth  mject  large  quantities  of  dust  and  debris 
high  into  the  atmosphere.  The  net  iniectton  m  a 
full-scale  nuclear  exchange  is  estimated  to  be  as 
much  as  several  hundred  teragrams  [Turco  et 
al.,  1983:  National  Research  Council  iNRC). 
1985:  Pittock  et  al..  1986:  Bacon.  Dunn,  and  Sar- 
ma.  19881  ■  The  presence  of  so  much  dust  in  the 
atmosphere  could  have  many  serious  conse¬ 
quences.  ranging  from  interference  with  military 
operations  to  changes  in  the  global  climate.  Sim¬ 
ple  models  have  been  developed  to  help  assess 
the  impact  of  nuclear  dust.  Yet.  because  our 
knowledge  of  the  lotting,  transport,  removal,  and 
optical  effects  of  large  atmospheric  dust  masses 
IS  limited,  these  models  rely  on  several  assump¬ 
tions  that  are  highly  uncertain. 

Opportunities  to  experimentally  study  the  physics 
of  nuclear  dust  clouds  ended  with  above  ground 
testing  in  1962.  Available  data  are  either  insuffi¬ 
cient  or  not  completely  analyzed.  One  alternative 
IS  to  identify  natural  events  that  display  some  of 
the  key  physical  mechanisms  believed  to  be  im¬ 
portant  in  nuclear  dust  clouds.  Two  natural  ana¬ 
logues  are  promising:  (1)  plumes  from  erupting 
volcanos  and  (2)  Saharan  dust  storms.  Volcanic 
plumes  have  been  considered  an  analogue  for 
high  altitude  dust  transport  [see.  for  example. 
NRC.  1985].  Such  plumes  may  also  model  con¬ 
densation  processes  (for  crater  ejecta)  that  oc¬ 
cur  at  upper  altitudes.  In  this  report,  wo  examine 
Saharan  dust  plumes,  and  evaluate  their  rele¬ 
vance  as  an  analogue  for  nuclear  dust  clouds. 

The  most  attractive  feature  of  Saharan  plumes 
is  the  amount  of  dust  involved.  A  single  storm 
usually  lasts  about  four  or  five  days,  and  typically 


transports  8  to  10  Tg  of  dust  westward  from  ’he 
Sahara  over  the  Atlantic  Ocean  [Carlson.  19'9| 
Assuming  a  dust  lofting  capacity  for  nuclear  sur¬ 
face  bursts  of  0.3  Tg/MT  [Rosenblatt.  1981].  a 
Saharan  dust  episode  equals  roughly  thirty  i  -MT 
surface  detonations  or  about  one-fifth  the  infec¬ 
tion  expected  from  a  full-scale  icounterforcei  at¬ 
tack  against  a  single  U  S.  Minuteman  missile 
field'  [Bacon.  Dunn,  and  Sarma,  1988). 

Since  Saharan  plumes  contain  such  a  large  quan¬ 
tity  of  lofted  material,  they  are  easily  monitored 
by  satellites.  Over  the  ocean,  the  highly  reflective 
plume  contrasts  well  against  the  water.  More¬ 
over,  aerosol  particles  of  true  Saharan  origin 
(which  are  usually  composed  of  quartz  grams  and 
clay)  and  locally  generated  aerosols  (i.e..  sea 
salt  particles)  are  readily  distinguished  at  sam¬ 
pling  stations  far  from  the  source  Thus,  key  as¬ 
pects  of  the  dust  plume  such  as  its  density  and 
its  panicle-size  distribution  (PSD)  can  be  mea¬ 
sured  as  a  function  of  downstream  distance.  To 
the  extent  that  a  Saharan  plume  resembles  a  nu¬ 
clear  dust  cloud,  we  can  learn  about  1 1 )  the  pos¬ 
sible  optical  and  climate  effects  of  nuclear  dust. 
(2)  radioactive  fallout  patterns.  (3)erosive  and 
mechanical  hazards  dust  presents  to  aircraft, 
and  (4)  sedimentation  rates. 

In  Sec.  2.  we  review  meteorological  and  geo¬ 
graphical  aspects,  long-range  transport  and  re¬ 
moval  mechanisms,  and  optical  effects.  The 
physical  mechanisms  involved  in  the  lofting  of 
wind-blown  dust  are  discussed  in  Sec.  3.  Saharan 
dust  plumes  as  analogues  for  nuclear  dust  clouds 
are  evaluated  in  Sec.  4,  and  our  conclusions  are 
summarized  in  Sec.  5. 


1  Assuming  185  targets  par  flald  and  a  scanarlo  involving  on#  500-KT  surfaca  burst  and  ona  500-KT  low-altituda 
alrburst  par  targal 


SECTION  2 


SAHARAN  DUST  PLUMES 


SOURCE  STRENGTH. 

The  Sahara  is  the  single  largest  source  of  wind¬ 
blown  atmospheric  dust  in  the  world.  It  is  esti¬ 
mated  that  the  Sahaia  generates  from  130  to  330 
Tg  per  year  [Prospero  and  Carlson,  1972:  Jae- 
nicke  and  Sch  1978:  Qanor  and  Mamane.  1982). 
This  IS  roughly  one-quarter  to  one-half  of  the  total 
annual  dust  generated  by  all  sources  worldwide, 
and  accounts  for  about  the  same  fraction  of  the 
total  dust  burden  found  in  the  atmosphere  at  any 
given  time  [Junge,  1979).  The  Sahara  probably 
contributes  almost  as  much  mass  to  the  global 
dust  budget  as  all  other  major  sources  combined 
te  g.,  the  Raiashan  Desert  of  India,  the  Gobi  Des¬ 
ert  of  Mongolia  and  northern  China,  and  the  des¬ 
erts  of  Mexico  and  the  southwestern  U.S.).  The 
Saharan  dust  mass  is  usually  added  to  the  atmo¬ 
sphere  if’  10  to  12  storms  per  year,  ranging  in 
Size  from  O.l  to  lOTg.  The  largest  storms  raise 
a  mass  comparable  to  surface  bursts  yielding 
tens  of  megaton.  The  source  region  for  Saharan 
dust  storms  can  cover  several  thousand  square 
kilometers,  roughly  equivalent  to  a  section  of  a 
missile  field. 

DISTRIBUTION  PATTERNS. 

We  are  concerned  primarily  with  the  westward 
transport  of  dust  from  the  Sahara  to  the  Atlantic 
Ocean.  Significant  amounts  are  also  carried  in 
other  directions.  Figure  1  shows  a  composite  of 
several  dust  transport  estimates.  Although  most 
of  the  lofted  mass  is  carried  westward,  approxi¬ 
mately  20  percent  (70  Tg/year)  is  transported 
eastward  toward  the  eastern  Mediterranean. 
Lesser  amounts  (in  unknown  quantities)  are  also 
carried  north  and  south,  as  evidenced  by  occa¬ 
sional  dust  falls  over  Europe  (Tullet,  1978:  Prodi 
and  Fea.  1979)  and  the  Gulf  of  Guinea  (Bertrand. 
Baudet.  and  Drochon,  1974). 

Most  of  the  dust  raised  from  the  Sahara  consists 
initially  of  particles  with  radii  greater  than  1  pm. 
These  particles  have  high  fall  velocities  and 
therefore  do  not  remain  lofted  for  very  long.  Fig¬ 
ure  1  shows  the  magnitude  of  this  prompt  deple¬ 
tion:  Of  the  approximately  330  Tg/year  raised 
over  the  Sahara,  only  about  one-quarter  (90  Tg/ 
year)  eventually  leaves  the  African  continent.  Of 


the  remainder,  about  40  percent  05  Tg  year'  s 
deposited  in  the  Atlantic  Ocean,  another  40  per¬ 
cent  crosses  the  Atlantic,  and  20  percent  tails 
over  the  Near  East.  Although  a  significant  fraction 
remains  in  the  atmosphere  and  is  transported 
several  thousand  kilometers  by  wind  systems, 
most  particles  fall  out  in  the  first  500  to  lOOO  km 
from  the  source.  The  depletion— mostly  m  micron 
or  larger  sized  particles— is  similar  to  that  antici¬ 
pated  in  nuclear  clouds. 

SEASONAL  DIFFERENCES. 

Westward  dust  transport  from  the  Sahara  is  per¬ 
ennial.  but  seasonal  differences  in  the  distribution 
occur  (caused  by  changes  in  major  wind  patterns 
across  Africa) .  Figure  2  shows  the  average  sum¬ 
mer  and  winter  locations  of  the  Saharan  plume 
from  satellite  observations  [Schutz.  1980).  The 
mean  streamlines  at  the  effective  steering  alti¬ 
tude  (about  3  km)  are  also  shown,  in  winter  iFig, 
2a) .  direct  transport  from  the  Sahara  to  the  Atlan¬ 
tic  is  blocked  by  the  mean  westerlies  over  north 
Africa.  The  plume  therefore  emerges  somewhat 
below  the  Sahara  over  the  Gulf  of  Guinea  and  ex¬ 
tends  southwestward  toward  Brazil  with  the  mean 
easterlies.  In  summer  (Fig.  2b).  these  easterlies 
shift  northward  and  direct  transport  from  the  Sa¬ 
hara  resumes.  The  summer  dust  is  then  earned 
almost  due  west.  This  seasonal  shift  is  therefore 
responsible  for  the  sharp  summer  maximum  (and 
winter  minimum)  in  dust  concentrations  noted  at 
Barbados  Island  [Carlson  and  Prospero,  1972). 

Although  these  seasonal  differences  are  mostly 
due  to  changes  in  weather  patterns,  there  are 
also  seasonal  differences  in  dust  cloud  composi¬ 
tion  (Delany  et  al..  1967:  Prospero  et  al..  1970). 
In  summer  (the  rainy  season  in  sub-Saharan  Afn 
ca).  grasses  inhibit  dust  scouring.  Similarly,  rain¬ 
fall  changes  the  soil  cohesiveness.  Consequent¬ 
ly.  the  seasonal  dependence  is  strong.  A  similar, 
but  more  pronounced  dependence  is  expected 
for  U.S.  and  Soviet  target  areas,  especially  in  re¬ 
gions  where  soil  freezing  occurs. 

PLUME  DEVELOPMENT. 

The  large  summertime  Saharan  dust  loads  are 
caused  by:  M)  the  presence  of  a  strong,  per- 
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a.  Winter. 


Figure  2.  Av/erage  Saharan  plume  positions  and  mean  700  mbar  (<\.  3  km)  streamline  fields  for  winter 
and  summer.  (Location  of  Sahara  indicated  by  hatching.) 
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sistent  easterly  let  stream  at  midtropo  spheric 
levels  and  (2)  intense  heating  of  the  desert  sur¬ 
face  by  solar  radiation.  As  air  masses  moving 
with  the  easterly  current  cross  the  desert,  they 
are  rapidly  warmed  from  below  by  contact  with 
the  heated  surface.  Widespread,  vigorous  dry 
convection  develops,  effectively  mixing  the  warm 
air  near  the  surface  to  higher  levels.  As  a  result, 
a  warm,  nearly  isentropic  turbulent  mixed  layer 
forms,  gradually  deepening  as  the  air  mass  pro¬ 
ceeds  westward.  As  this  deepening  continues, 
some  of  the  horizontal  momentum  contained  in 
the  easterly  jet  is  carried  to  lower  levels  through 
the  neutrally  stable  mixed  layer.  Surface  wind 
speeds  increase  and  dust  lofting  ensues.  The  re¬ 
sult  is  a  Saharan  dust  storm.  (The  precise  condi¬ 
tions  necessary  for  the  initiation  of  dust  lofting 
from  the  surface  are  discussed  in  Sec. 3.) 

The  same  convective  thermals  that  trigger  a  dust 
storm  also  mix  the  dust  vertically.  This  vertical 
transport  continues  until  the  air  mass  reaches  the 
Atlantic.  The  mixed  layer  is  often  completely  dust 
filled  and  can  extend  to  6  km  in  altitude  [Carlson 
and  Prospero.  1972).  Potential  temperatures 
within  the  layer  are  quite  high  (-316  or  317  K  in 
July),  and  nearly  constant  with  height.  During 
summer,  this  region  is  usually  dominated  by  a 
strong  high-pressure  system,  and  the  air  above 
the  layer  can  be  warmer.  Consequently,  a  weak 
capping  inversion  forms,  effectively  trapping  the 
dust  below  it  and  inhibiting  further  deepening  of 
the  mixed  layer. 

Saharan  plume  dust  concentrations  over  land  are 
roughly  700  ug/m^  [Westphal.  Toon,  and  Carl¬ 
son.  1907].  Assuming  a  constant  mass  concen¬ 
tration  to  4  km.  the  total  column  density  is  2800 
kg/ km2— roughly  250  times  the  average  Northern 
Hemisphere  background  dust  burden  [Junge. 
1979).  The  source  areas— chiefly  dry  lake  and 
stream  beds— although  relatively  small  (on  the 
order  of  10^  km2),  are  distributed  over  a  very 
large  region  [Schutz.1980] .  The  dust  from  these 
multiple  sources  merges  into  a  single  mesoscale 
plume.  A  schematic  of  the  early  dust  plume  de¬ 
velopment  and  its  evolution  over  the  Atlantic  is 
shown  in  Fig. 3. 

SAHARAN  AIR  LAYER. 

Once  over  the  ocean,  the  Saharan  dust  plume 
encounters  the  northeast  trade  winds.  These 
winds  are  almost  always  associated  with  a  strong 


temperature  inversion  based  between  0  5  and 
1.5  I'm  altitude  (Riehl.  1979).  The  manno  air 
trapped  below  this  inversion  is  initially  dust  ''ee 
and  is  much  cooler  (and  denser)  than  the  Saha¬ 
ran  air.  Consequently,  as  the  dusty  air  contained 
in  the  deep  mixed  layer  emerges  from  the  conti¬ 
nent.  it  is  forced  above  the  trade  wind  inversion 
(Fig.  3) .  The  elevated  air  mass  is  called  the  Saha¬ 
ran  air  layer  (SAL). 

The  temperature  sounding  in  Fig,  4  shows  that 
the  SAL  is  distinguished  by  a  thick  elevated  layer 
of  nearly  isentropic  air  in  the  lower  troposphere. 
The  SAL  is  well  mixed— a  result  of  its  convective 
origin  over  the  Sahara,  which  accounts  for  the 
near-uniform  vertical  distribution  of  fine  ( submi¬ 
cron  radius)  dust  particles  [Prospero  and  Carl¬ 
son.  1972).  The  vertical  extent  of  the  mixing  is 
limited,  however,  by  the  presence  of  two  temper¬ 
ature  inversions:  the  trade  wind  inversion,  which 
marks  the  lower  boundary,  and  the  capping  inver¬ 
sion.  based  around  6km.  which  defines  the  upper 
boundary.  Since  these  inversion  layers  are  highly 
statically  stable,  vertical  mixing  through  them  is 
inhibited.  Fine  dust  is  therefore  trapped  between 
the  inversions  and  “fumigates"  from  the  SAL  into 
the  underlying  marine  layer  at  a  very  slow  rate. 
Consequently,  dust  concentrations  within  the 
marine  layer  are  usually  two  to  four  times  lower 
than  they  are  in  the  SAL  [Prospero  and  Carlson, 
1972:  Diaz  et  al..  1976).  Thus  the  SAL  channels 
the  dust  in  a  well-defined  layer  across  the  Atlan¬ 
tic. 


The  SAL  is  a  prime  example  of  how  interaction 
with  a  weather  system  can  influence  long-range 
dust  transport.  Nuclear  dust  clouds  originating  m 
mid-latitudes  will  be  similarly  influenced.  For  ex¬ 
ample.  dust  masses  passing  over  the  oceans 
could  be  isolated  from  near-surface  turbulent 
mixing  and  removal  by  low-altitude  inversions, 
especially  in  summer.  On  the  other  hand, 
large-scale  organized  precipitation  systems  are 
more  prevalent  in  mid-latitudes  than  in  the  trop¬ 
ics.  It  is  therefore  likely  that  a  nuclear  cloud  will 
encounter  one  of  these  systems  within  a  few 
thousand  kilometers  of  its  source.  Prompt  re¬ 
moval  of  the  dust  by  precipitation  scavenging  (at 
least  in  the  lower  troposphere)  would  then  ensue. 
This  mechanism  will  be  enhanced  if  the  cloud 
contains  large  concentrations  of  submicron-sized 
clay  particles  (especially  kaolinite.  illite.  and 
montmorillonite) .  Such  particles  serve  as  effi¬ 
cient  ice  nuclei  at  temperatures  as  high  as  -3’C 
[Mason.  1971). 
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Figure  3.  Schematic  of  ma|or  processes  involved  in  lofting,  transport,  and  lemoval  ol  Saharan  dust  liom  aimosptieie 
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Figure  4.  Typical  summer  temperature  profile  over  tropical  north  Atlantic  near  Aliican  coast 


MASS  DEPLETION  RATES. 

In  the  absence  of  in  situ  depletion  measure¬ 
ments.  we  have  analyzed  the  vertical  dust  density 
protiles  predicted  by  two  simulations  of  the  Saha¬ 
ran  dust  plume  [Lee.  1983;  Schutz.  1980).  By 
integrating  the  profiles  vertically,  the  change  in 
column  density  M  (and  hence  the  plume  deple¬ 
tion  rate  R)  was  determined  as  a  function  of  dis¬ 
tance  downstream  d.  Figure  5  shows  M(d) 
expressed  as  a  percent  of  the  initial  column  den¬ 
sity  Mq. 

Clearly,  the  depletion  rate  decays  with  distance 
from  the  source.  Based  on  Fig.  5.  we  estimate 
that  roughly  one-half  the  plume  mass  is  depleted 
in  the  first  1000  km  and  that  about  70  percent  falls 
out  within  200  km.  The  column  density  is  given 
by 


Mamane  (I982j  and  the  optical  depth  piume 
mass  measurements  of  Carlson  (1979]. 

From  Eqs.(1)  and  (2).  the  plume  depletion  rate 
following  a  downstream  trajectory  is 


(A  +  Bd)  2 


'3) 


where  U  is  the  horizontal  transport  speed.  As¬ 
suming  Mq  =  3000  kg/km2  and  U  =  20  m/s  (typical 
values  for  the  Saharan  plume) .  we  find  R  decays 
exponentially  from  an  initial  value  Rg  =  236  kg- 
km2/h  as  shown  in  Fig. 6.  Far  from  the  source  'd 
7000  km),  the  depletion  is  about  two  or- 
ders-of-magnitude  smaller  than  Rq  (R  SI  2kg 
km2/h):  this  rate  corresponds  to  a  surface 
deposition  of  only  a  few  microns  per  year. 


M(d)  =  Mo  0  ,  (1) 

where  the  characteristic  “e-folding"  distance  D 
is  itself  a  function  of  d.  We  find  that  to  a  very  good 
approximation 

D(d)  =  A  +  Bd  ,  (2) 

where  A  =  91 5  km  and  B  =  0.34  km/km  in  the  Sa¬ 
haran  plume.  Thus  O  increases  linearly  down¬ 
stream.  This  is  somewhat  different  than  the 
predictions  of  Lee  (19831  and  Schutz  (19801.  In 
those  models,  the  plume  was  assumed  to  origi¬ 
nate  either  1400  km  inland  over  the  Sahara 
[Lee.  1983).  or  on  the  coast  (Schutz.  1980). 
Thus,  the  meteorological  conditions  encountered 
by  the  simulated  plumes  were  quite  different  dur¬ 
ing  the  early  stages  of  their  evolution.  Moreover. 
Schutz  and  Lee  assumed  different  initial  dust 
concentrations.  PSOs.  and  horizontal  transport 
speeds  in  their  plumes,  and  treated  vertical  mix¬ 
ing  and  surface  removal  differently.  In  particular. 
Lee  (1983)  assumed  the  surface  deposition  ve¬ 
locity  to  be  independent  of  particle  radius.  We 
therefore  estimate  a  somewhat  faster  depletion 
than  Lee  (1983)  in  the  first  4(XX3  km  and  slower 
at  greater  ranges  (see  Fig. 5) .  This  rate  is  consis¬ 
tent  with  the  mass  flux  estimates  of  Qanor  and 


PARTICLE-SIZE  DISTRIBUTION. 

The  radiative,  optical,  and  mechanical  hazard 
properties  of  any  dust  plume  depend  not  only  on 
the  total  mass  contained  within  the  plume,  but 
also  the  PSD.  Windblown  or  scoured  dust  usually 
comprise  a  wide  range  of  particle  sizes. 

Oust  distributions  representative  of  near-surface 
PSDs  in  the  source  region  (Libya)  and  at  three 
points  downstream  are  shown  in  Fig. 7.  The  cu¬ 
mulative  mass  per  unit  volume  m  of  all  particles 
of  radii  less  than  r  Is  given  as  function  of  particle 
size. 2  Figure  7  shows  that  (1 )  the  total  mass  con¬ 
centration  decays  very  rapidly  downstream.  (2) 
with  distance  from  the  source,  the  distribution 
shifts  to  smaller  particles,  and  (3)  the  depletion 
rate  depends  strongly  on  particle  size.  This  is  es¬ 
pecially  obvious  when  we  compare  the  relative 
change  in  dm/d  (log  r)  for  radii  loss  than  1  jim  with 
that  for  larger  radii.  Most  of  the  mass  depletion 
occurs  wiren  larger  particles  (r  >  1  4m)  fall  out. 
The  submicron  concentration  remains  fairly  con¬ 
stant  (Table  1);  these  particles  are  transported 
long  distances— even  though  the  plume  top  is 
fairly  low.  Nuclear  clouds  rising  to  the  upper  tro¬ 
posphere  would  travel  even  farther  and.  if  they 
penetrate  into  the  stratosphere,  could  remain 
aloft  for  years.  (Volcanic  plumes  provide  a  good 
analogue  for  such  stratospheric  transport.) 


2.  Th«  mass  dtstrHjutlon  dm/d|log  rj  Is  proportional  to  tha  numoar  dlstiibutton  dN/dr  by  dN/dr*  (3/4-tir‘poiniO) 
[dm/d|logr|| .  wttara  pg  Is  tha  particia  dansity. 
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Normalized  column  density  in  Saharan  plume  as  lunction  of  distance  downstieam 
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Figure  6.  Mass  depletion  rate  in  Saharan  plume  as  lunction  of  distance  downstieam 


Representative  particle-size  distiibutions  in  Saliaian  plume 


Table  1.  Mass  partitioning  in  Saharan  dust  observations. 


Observation 

Site 

Total 

Mass 

Concentration 

(ug/m3) 

Percent  by 

Mass  with 

Radii  Less 

Than  1  p.m 

Mass  Concentra¬ 
tion  of  Sub¬ 
micron  Portion 

(4g/m3) 

Median 

Radius 

lumi 

Central  Sahara 

W2800 

~  1 

~28 

~  38.6 

Cape  Verde  Islands 

215 

6 

13 

8.9 

Central  Atlantic 

75 

11 

3.3 

4.2 

Western  Atlantic 

5 

44 

2.2 

2.8 

Saharan  plumes  are  wind  lofted  (scoured)  and  rain  systems  are  more  frequent  and  wet  scav- 

therefore  consist  mostly  of  smaller  sired  (r  <  1  enging  in  the  troposphere  is  likely, 

nm)  particles— at  least  at  later  times.  The  scour¬ 
ing  mechanism  for  nuclear  clouds  is  much  stron-  Coagulation  has  a  noticeable  impact  on  the  evo- 

ger  (although  transient)  and  initially  a  wider  lotion  of  the  Saharan  PSD  only  when  turbulent 

spectrum  of  particle  sizes  is  lofted.  However,  the  mixing  and.  especially,  sedimentation  are  not 

largest  particles  fall  out  promptly.  At  later  times.  dominant  (Lee.  1983).  Rapid  depletion  by  sedi- 

nuclear  cloud  PSOs  may  therefore  be  similar  to  mentation  and  turbulent  mixing  lessens  the  devei- 

the  Saharan  distributions.  opment  of  large  particles  by  coagulation  (see  Fig. 

8).  For  Saharan  plumes,  coagulation  is  usually 
REMOVAL  MECHANISMS.  not  a  significant  factor.  There  have  been  mea¬ 

surements  made  over  the  Atlantic,  however,  that 
A  dust  particle  can  be  removed  by  sedimenta-  show  particles  too  large  to  have  survived  deple¬ 
tion,  mixing  to  the  surface  by  turbulence,  scav-  tion  by  sedimentation  and  turbulent  mixing.  Par- 

enging  by  precipitation  (wet  removal),  or  by  tides  in  excess  of  1(X)-nm  radius  have  been 

combining  with  one  or  more  particles  to  form  a  sampled  in  the  marine  layer  over  the  Canary  is- 

larger  aggregate  (coagulation).  At  least  two  of  lands  [Goude-Gaussen  et  al. .  1987] .  Coagulation 

these  mechanisms— sedimentation  and  turbulent  obviously  is  responsible.  Observations  of  anoma- 

mixing— play  key  roles  in  Saharan  plumes.  Sedi-  lously  large  particles  have  also  been  made  over 

mentation,  which  depends  on  the  terminal  vetoc-  the  Cape  Verde  Islands  (QIaccum  and  Prospero. 

ity  of  a  particle  and  hence  on  its  size,  is  1980]  and  the  central  Atlantic  (Carder  et  al.. 

responsible  for  the  rapid  removal  of  large  (r  >  1  1986). 

pm)  particles.  Turbulent  mixing  is  effective  for 

particles  so  small  that  their  sedimentation  is  sig-  Clearly,  some  mechanism  other  than  direct 

nificantly  retarded  by  viscous  forces  and  is  re-  transport  from  the  Sahara  is  responsibie  for  such 

sponsibie  for  the  gradual  removal  of  submicron-  “ uUragiant "  particles.  One  possibility  is  local 

radius  particles.  “cloud  processing,"  which  may  occur  when 

small,  initially  uncoagulatod  particles  from  the 
Wet  removal  apparently  has  little  effect  on  Saha-  SAL  pass  through  underlying  stratus  decks 

ran  plumes,  since  the  transport  occurs  mostly  in  (Westphal.  Toon,  and  Carlson.  1987).  Such  a 

the  dry  air  between  successive  westward  propa-  moisture-aided  coagulation  mechanism  would  be 

dating  wave  disturbances  (Carlson  and  Prosoero.  important  when  dust  clouds  encounter  nonpre- 

1972],  This  is  quite  different  from  the  situation  cipitating  clouds.  But  this  is  speculative:  the  pres- 

expected  for  mid-latitude  nuclear  clouds,  where  ence  of  ultragiant  particles  far  downstream  re- 


Figure  8.  Effect  of  various  dust  removal  mechaitisins  on  Saharan  dust  pai tide  size  distnbuiioii  atlei  48  h  simulation  time 


mains  unexplained.  The  possibility  that  large  par¬ 
ticles  could  exist  far  from  the  source  is  relevant 
for  nuclear  clouds  since  these  particles  lespe- 
cially  for  r  >  20  umi  could  cause  the  most  me¬ 
chanical  damage  to  flight  vehicles.  Moreover,  in 
the  mid-latitudes,  the  opportunity  for  formation 
of  larger  particles  by  cloud  processing  is  probably 
much  greater  than  in  the  Saharan  dust  systems. 

RADIATIVE  EFFECTS. 

Dust  is  an  obscurant.  A  beam  or  signal  through 
a  cloud  of  dust  is  attenuated  by  absorption  and 
scatteiing.  The  amount  of  this  attenuation  is  giv¬ 
en  by  Beer’s  law,  which  expresses  the  ratio  of 
the  beam  irradiance  lx  at  wavelength  ,\  emerging 
from  the  cloud  to  the  incident  irradiance  l\o. 

('x/'v)"*' 

Here  t\  is  the  spectral  optical  depth.  For  a  verti¬ 
cally  oriented  beam,  it  is  defined  as 


r»h 


where  z  is  the  vertical  coordinate,  h  is  the  cloud 
thickness,  and  kx  and  kx.j  are  the  volume  spec¬ 
tral  absorption  and  scattering  coefficients,  re¬ 
spectively. 

In  Saharan  plumes,  the  particles  are  sufficiently 
large  compared  to  solar  wavelengths  that  the  op¬ 
tical  depth  is  linearly  related  to  the  column  mass 
M  by 

T  =  M  A*  (6) 

,v 

where  iy  is  constant  equal  to  3750  kg/km2  near 
Africa  and  roughly  1500  kg/km2  over  the  western 
Atlantic  (Carlson  and  Caveriy,  1977].  This  formu¬ 
la  was  used  to  derive  the  satellite-observed  col¬ 
umn  densities  shown  in  Fig.  5. 

Figure  9  maps  tx  at  a  solar  wavelength  (x  =  0.66 
nm)  determined  from  satellite  observations  dur¬ 
ing  a  particularly  severe  Saharan  dust  episode. 
Note  the  elongated  pattern  of  the  isooleths.  The 
largest  values  (tx  >  3.5)  are  located  near  the 
Mauritanian  coast  where  the  olume  first  emerges. 
As  the  dust  is  transported  westward  and  spreads 


laterally,  measured  optical  depths  lessen,  con¬ 
sistent  with  disoersion  and  deoietmn  of  mo 
plume.  Values  of  tx  m  excess  of  2.0  co'^er  on 
area  of  about  3  x  10®  km2.  These  optical  deoths 
are  very  high— oackground  values  m  relatively 
dust-free  regions  seldom  exceed  0.25— and  solar 
irradiance  at  the  surface  is  reduced  by  36  per¬ 
cent.  Such  obscurations  are  consistent  with  sur¬ 
face  dust  concentrations  of  200  ug;m3  measured 
at  Sal  Island  (in  the  Cape  Verde  islands) .  At  these 
concentrations,  visibility  is  reduced  to  the  point 
that  it  is  impossible  to  distinguish  cloud  from  sky 
[Carlson,  1979).  The  outbreak  shown  m  Fig. 9 
was  exceptionally  massive;  a  more  typica'  value 
of  TX  in  Saharan  plumes  is  1.5  [Carlson.  1979] . 
A  nuclear  cloud,  however,  would  be  at  least  twice 
as  thick  as  the  Saharan  plume  at  these  early 
times.  The  optical  depth  would  also  be  twice  as 
large.  The  large  outbreak  shown  in  Fig.  9  is  there¬ 
fore  more  representative  of  the  optical  depths  ex¬ 
pected  in  a  nuclear  dust  cloud. 

The  large  optical  depths  associated  with  some 
Saharan  dust  storms  imply  significant  modifica¬ 
tion  of  local  weather,  especially  near  the  source 
region.  Even  several  hundred  kilometers  south  of 
the  Sahara  (Nigeria)  visibilities  during  winter  dust 
eoisodes  can  drop  to  less  than  0.5  km  at  the  sur¬ 
face  [Brinkman  and  McGregor.  19831  During 
these  episodes,  visibility  tends  to  be  low  during 
the  day  and  improves  at  night.  Thermal  convec¬ 
tion  from  the  heated  ground  helps  fumigate  lofted 
dust  to  the  surface.  Reductions  in  surface  tem¬ 
peratures  are  associated  with  the  low  visibilities. 
Daily  maximum  temperature  changes  up  to  6'’C 
have  been  reported  [Brinkman  and  McGregor. 
1983). 

In  Saharan  plumes  absorption  plays  a  minor  role 
compared  to  scattering.  The  magnitude  of  the 
absorption  is  determined  by  the  molecular  struc¬ 
ture  and  composition  of  the  dust  particles,  and 
is  measured  by  the  imaginary  part  of  the  index 
of  refraction  n|.  This  dimensioniess  quantity  'S  re- 
iated  to  the  volume  absorption  coefficient  f'x  a  by 


For  Saharan  dust,  nj  is  very  small;  Measured  val¬ 
ues  at  solar  wavelengths  range  from  about  O.OOi 
to  0.02.  with  an  average  around  0.005  [Carlson 
and  Benjamin.  1980).  Such  values  are  typical  of 
silicates.  By  comparison,  smoko,  which  is  a 
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much  more  efficient  adsorber  of  solar  radiation, 
has  an  imaginary  index  of  refraction  around  0. 1 
or  larger  [NRC.  1985] .  The  low  absorotivity  of  Sa¬ 
haran  dust  is  due  m  part  to  its  very  low  organic 
content.  Saharan  aerosols  contain  less  than  2 
percent  mass  of  carbonaceous  material  (Jae- 
nicke  and  Schutz.  1978] .  This  is  much  lower  than 
the  organic  fraction  found  in  most  other  dust 
aerosols  around  the  world  (supporting  the  view 
that  the  Sahara  is  ‘organically  dead.")  Nuclear 
dust  clouds  are  likely  to  be  more  absorptive  than 
Saharan  plumes,  since  many  U.S.  and  USSR  silo 
fields  are  located  on  thick,  organically  rich  praine 
or  steppe  soils  imollisols  or  chernozem) .  Stron¬ 
ger  absorption  in  these  clouds  would  increase  at¬ 
mospheric  warming  in  the  upper  troposphere  and 
could  induce  a  “self-iofting "  mechanism  similar 
to  that  expected  for  smoke  clouds.  This  has  not 
been  detected  m  the  largely  silicate  Saharan 
plumes. 


Scattering  in  dust  Clouds  is  a  strong  function  of 
particle  size.  In  the  Saharan  plume,  the  effective 
scattering  radius  r,  is  usually  quite  large  com¬ 


pared  to  solar  wavelengths  O.Siimi. 

example,  over  the  eastern  Atlantic,  r,  2  ’n  5 
pm  '.Carlson  and  Caverly,  1977] .  As  a  result,  'Ue 
scattering  coefficient  k)^  ,  is  large,  and  most  of 
the  incident  solar  radiation  is  reflected  or  back- 
scattered.  Further  downstream,  however,  the 
large  particles  responsible  for  reflection  are  de¬ 
pleted  and  rg  ft  Xpagt,.  Scattered  radiation  is  then 
directed  primarily  m  the  forward  direction  Mie 
scattering) :  k\  ,  therefore  decreases  and  the  op¬ 
tical  depth  due  to  scattering  s  reduced. 

In  a  nuclear  dust  cloud,  scattering  depenris  "n 
the  characteristic  particle  size.  For  a  doud  v,ith 
an  initial  PSD  similar  to  a  Saharan  distribution. 
expect  that  k\  ,  at  solar  wavelengths  would  de¬ 
crease  downstream  as  large  particles  are  de¬ 
pleted  and  the  scattering  passes  from  the 
geometric  to  the  Mie  regime,  it  is  possible,  how¬ 
ever.  that  nuclear  clouds  would  have  a  greater 
initial  mass  fraction  of  submicron-size  particles 
due  to  condensation  of  vaporized  soils.  Mie  scat¬ 
tering  would  then  dominate,  and  the  nuclear 
cloud  would  therefore  be  less  reflective  than  the 
Saharan  plume. 
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SECTION  3 


LOFTING  MECHANISMS 


VERTICAL  MASS  FLUX. 

Winds  at  the  earth's  surface  exert  shear  stresses 
Oh  the  ground.  If  the  stress  is  large  enough,  ma¬ 
terial  IS  raised  from  the  surface  and  transported 
by  the  airstream  (saltation).  The  vertical  mass 
flux  F,  is  proportional  to  the  wind  speed  as 

Fz  =  CU."  .  (8) 

where  C  and  n  are  empirical  constants  and  U*  is 
the  friction  velocity  for  a  turbulent  shear  layer 
[Gillette.  Blifford.  and  Fenster.  1972]  computed 
from: 

U,  =  k^U2-Uij  inj^Z,  Zaj  .  (9) 

Here  k  SI  0.4  is  the  Von  Karmann  constant  and 
Ui  and  Uj  are  wind  speeds  measured  at  heights 
Zi  and  Z2.  respectively. 

The  experimentally  measured  dependence  of  F^ 
on  U*  for  three  desert  soil  types— sand.  loam, 
and  clay  is  summarized  in  Fig.  10.  The  vertical 
mass  flux  increases  nonlinearly  with  increasing  U* 
and  depends  on  soil  type.  A  least-squares  fit  by 
Westphal.  Toon,  and  Carlson  [1987]  suggests 
that 

Fz  a  2x10-'Ou3  kg/m2/s  ,  (IO) 

where  U*  is  in  centimeters  per  second.  Thus  n 
=  3  in  Eq.(8).  The  vertical  mass  flux  of  particles 
from  the  surface  is  therefore  proportional  to  the 
flux  of  kinetic  energy  delivered  to  the  surface  [Gil¬ 
lette.  1980].  However,  shock  tube  experiments 
[Hartenbaum.  1971]  and  theoretical  work  [Mi- 
rels.  1984]  show  that  pa  1.1  to  1.2  at  wind 
speeds  more  typical  of  a  nuclear  environment  [U* 
=  0(102)  m/sj.  Both  Hartenbaum's  [1971]  re¬ 
sults  and  Mirels's  [1984]  are  for  highly  idealized 
conditions  (e  g.,  flat,  noncohesive  soil  surfaces) 
rattier  than  natural  soils.  The  greater  scouring 
rates  for  Saharan  dust  may  also  result  from  the 
steady  winds,  whereas  in  a  nuclear  pulse,  peak 
velocities  are  large  but  transient. 

The  measurements  in  Fig.  10  also  show  a  thresh¬ 
old  friction  velocity  U-jh  below  which  no  erosion 


occurs.  Moreover,  the  value  of  U->h  deoenris 
strongly  upon  the  soil  type,  as  indicated  by  the 
vertical  dashed  lines  in  Fig.  10. 

THRESHOLD  VELOCITY. 

Soil  types  and  conditions  vary  widely  deoeoding 
on  location,  time  of  year,  and.  in  deserts,  the 
diurnal  cycle.  The  threshold  friction  velocity 
quired  to  initiate  dust  lofting  also  varies,  although 
the  dependence  of  U-th  soil  type  is  not  imme¬ 
diately  obvious.  For  example,  clay  soils  ‘which 
contain  a  large  number  of  very  tine  submicron 
particles)  have  higher  threshold  velocities  than 
sand.  The  distinction  Is  important  since  U-th  de¬ 
termines  whether  dust  storms  develop. 

Values  of  U*th  for  a  variety  of  soil  types  and  con¬ 
ditions  in  the  Mojave  Desert  have  been  exper¬ 
imentally  determined  [Gillette  et  al. .  1980],  Both 
undisturbed  and  disturbed  soils  were  tested.  The 
results  are  presented  in  Table  2.  Disturbed  soils 
generally  have  lower  threshold  velocities  than  un¬ 
disturbed  soils:  the  presence  of  clay  or  salts 
tends  to  increase  U-th'  especially  if  a  crust  is 
present.  Even  a  thin  (0.5-cm  thick)  and  weak 
(modulus  of  rupture  0.7  bar)  crust  (e.g.,  typi¬ 
cal  of  undisturbed  dry  lake  beds)  can  protect  the 
underlying  soil  from  all  but  the  strongest  winds 
[Gillette  et  al..  1982].  The  strength  of  a  crust  is 
highly  dependent  on  the  amount  of  mica-free  clay 
present,  although  exchangeable  sodium  and  cal¬ 
cium  carbonate  also  increase  cohesiveness. 
Scouring  thresholds  are  also  increased  for  freeze 
conditions  and  snow  coverings,  water  or  salt  (es¬ 
pecially  in  crystalline  form)  in  the  soil,  aggrega¬ 
tion  or  “clumping"  of  particles,  a  high  clay 
content,  and  obstacles  which  increase  surface 
roughness  (e.g.,  peb-  bles.  boulders,  vegeta¬ 
tion)  . 

Modifying  the  threshold  velocity  to  account  for 
the  nuclear  environment  is  clearly  necessary.  In¬ 
tense  thermal  radiation  can  fuse  soil  particles  and 
form  a  thin  but  strong  glassy  crust.  On  the  other 
hand,  thermal  radiation  also  vaporizes  soil  water, 
leading  to  a  decrease  in  cohesiveness  and  m 
U-th-  ih  some  cases,  the  rapid  vaporization  is  ex¬ 
plosive.  and  can  disrupt  a  preexisting  crust. 
Shock  loading  at  the  surface  can  also  break  soil 
crusts  (spalling)  and  produce  large  aggregated 
soil  clumps.  Such  clumps  have  much  larger  aero- 
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Figure  10.  Dependence  of  vertical  dust  mass  flux  tiom  surlace  on  Inciion  velocUy 


Table  2. Threshold  velocities  for  various  soil  types. 


Soil  Type 

Approximate 

Threshold  Velocity  Range 

(m/si 

Disturbed  soils  with  <  50%  clay 

0.20-0.40 

content  and  sparse  pebble  cover 

(<  20%  surface  area) 

Tilled  bare  sand 

0.25-0.40 

Disturbed  pebbly  soils 

0.40-0.65 

Bare  clay  with  naturally  disaggregated 

0.45-0.70 

particles 

Disturbed  soils  with  high  salt  or  clay  content 

0.40-1.50 

Undisturbed  sandy  soils 

1.40-2.00 

Undisturbed  soils  with  high  clay  content 

>  1.5 

and  surface  crusts 

Salt-crusted  soils 

>  1,5 

Soils  covered  by  coarse  (>  5  cm)  pebbles 

>  1.8 

Source:  Gillette  [1980). 


dynamic  cross  sections  than  individual  soil  grains 
and  are  therefore  preferentially  lofted  from  the 
surface  (assuming  large  enough  U*).  However, 
they  are  massive,  do  not  remain  lofted  for  long. 


and  break  into  smaller  pieces  at  the  surface  [Pat¬ 
terson  and  Gillette.  1977) ,  The  soils  can  then  be 
scoured  and  carried  to  high  altitude  with  the  rising 
fireball. 
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SECTION  4 


NUCLEAR  DUST  ANALOGUE 


The  Saharan  dust  plume  transports  large  dust 
masses  thousands  of  kilometers,  although  m  a 
cloud  of  limited  vertical  extent.  Figure  t1  com¬ 
pares  predicted  maximum  dust  concentrations  at 
ten  hours  for  two  nuclear  exchange  scenarios-^ 
with  that  for  a  natural  Saharan  plume.  The  Saha¬ 
ran  plume  IS  clearly  a  lower  tropospheric  phe¬ 
nomenon.  while  nuclear  clouds  involve  the  entire 
troposphere  and  part  of  the  stratosphere.  Never¬ 
theless.  the  dust  concentrations  within  nuclear 
and  Saharan  dust  clouds  are  similar  at  low  alti¬ 
tudes  (below  3  km). 

Other  aspects  of  the  Saharan  plume  are  com¬ 
pared  with  nuclear  clouds  in  Table  3.  Clearly,  the 
analogue  is  not  perfect;  Soils  in  the  source  re¬ 
gions  are  dissimilar,  nuclear  sweepup  is  more  en¬ 
ergetic  (and  may  involve  blast  effects),  and 
removal  mechanisms  are  different  (especially  in 
middle  latitudes) .  There  are  similarities  however, 
especially  in  the  long-range  transport  and 
long-time  cloud  depletion. 


the  soil  type  and  its  condition  le.g..  whether  or 
not  it  is  crusted,  the  amount  of  day  and  salt,  the 
presence  of  vegetation).  Moreover,  the  total 
mass  lofted  and  fraction  transported  long  dis¬ 
tances  depends  on  the  soil  size  distribution. 
These  factors  are  clearly  important  for  the  weak 
Saharan  forces,  but  may  also  influence  the  char¬ 
acteristics  of  lofted  nuclear  dust  clouds. 

Silo  fields,  located  primarily  in  regions  dominated 
by  thick,  organically  rich  soil  deposits  le.g..  the 
Great  Plains) .  experience  winter  freezes  and  sea¬ 
sonal  changes  in  vegetation.  The  susceptibility  of 
such  target  area  soils  to  sweep-up  is  clearly  quite 
different  from  that  of  arid  (Saharan)  soils.  Never¬ 
theless.  large  dust  clouds  are  swept  up  by  sur¬ 
face  winds  in  both  cases:  and  far  from  the 
source,  the  depletion  can  be  similar. 

LONG-RANGE  TRANSPORT.  THE  DEVELOP¬ 
MENT  OF  THE 


LOFTING  MECHANISMS. 

Low-altitude  nuclear  explosions  loft  dust  into  the 
atmosphere  by  one  of  two  mechanisms:  (I) 
sweepup  by  the  strong,  blast-induced  winds  and 
(2)  cratering.  For  detonations  with  a  scaled 
height  of  burst  (SHOB)  less  than  about  70  ft/ 
MT'  d  tjoth  of  these  mechanisms  contribute  to 
the  lofted  dust  mass.  For  greater  SHOBs,  only 
sweepup  is  a  factor  (NRC,  1985). 

The  physics  of  sweepup  are  poorly  understood. 
Generally,  experimental  studies  have  been  con¬ 
ducted  under  highly  idealized  laboratory  situa¬ 
tions  (e.g..  Hartenbaum,  1971).  Natural  soils  and 
conditions  have  not  been  extensively  studied  un¬ 
der  nuclear  sweepup  circumstances  (thermal  ir¬ 
radiation  of  the  surface,  high  blast  overpres¬ 
sures.  high-velocity  transient  winds,  etc.).  Natu¬ 
ral  dust  lotting  processes  provide  a  simulation, 
if  only  partial,  of  nuclear  sweepup. 

Saharan  dust  clouds  vary  with  the  source  region 
and  time  of  year.  The  rate  at  which  dust  is  lofted 
or  scoured  from  the  surface  strongly  depends  on 


Saharan  plume  requires  specific  atmospheric 
conditions.  These  conditions  also  make  it  possi¬ 
ble  for  the  dust  to  be  carried  long  distances. 
Without  the  channeling  influence  of  the  Saharan 
air  layer  over  the  tropical  Atlantic  and  in  the  ab¬ 
sence  of  wet  removal,  it  is  doubtful  that  much 
dust  would  survive  in  the  atmosphere  as  long  as 
it  does.  Nuclear  dust  clouds  are  likely  to  be  sub¬ 
ject  to  a  full  range  of  meteorological  conditions. 
In  midlatitudes,  dust  depletion  rates  are  higher 
principally  due  to  organized  precipitation  sys¬ 
tems.  Thus,  the  Saharan  plume  probably  repre¬ 
sents  an  upper  limit  to  dust  cloud  transport  and 
longevity  in  the  lower  troposphere. 

DEPLETION  RATES. 

Although  there  are  obvious  differences,  far  from 
the  source  region  (or  target  area),  mass  deple¬ 
tion  for  Saharan  and  nuclear  dust  clouds  should 
be  similar.  For  the  two  cases,  the  constants  m 
Eq.(1)  are  different .  Nuclear  burst  lofting  mecha¬ 
nisms  (scouring  and  fireball  upsweep)  are  much 
stronger  than  in  the  Saharan  case,  and  the  initial 
lofted  mass  is  dominated  by  more  massive  par- 


3.  Scsnarto  A  assumas  a  U  S.  laydown  on  ttim  Soviat  Union  Involving  1500  MT  total  ylafd  with  500  MT  In  furfaca 
bursts.  In  scenario  B.  a  Soviet  counterlorce  attack  on  the  U.S.  Is  assumed  Involving  several  thousand  mega¬ 
tons  of  surface  and  near-surface  bursts,  including  several  large  (20MT)  detonations  (see  Yoon  at  ai.  19861 
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Figure  1 1 .  Gompanson  ot  simulated  vertical  dust  concentration  proliles  in  Saharan  plume  and  nuclear  dust  cloud 
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tides.  The  initial  depletion  is  there  fore  larger. 
Moreover,  the  largest  particles  are  rapidly  de¬ 
pleted  so  that  the  column  density  M  decays  more 
rapidly.  Thus,  the  e-folding  distance  D  is  smaller. 
Finally,  far  downstream,  where  depletion  is  domi¬ 
nated  by  submicron-size  particles.  M  is  likely  to 
be  larger  due  to  the  greater  Initial  submicron-sire 
concentrations  in  the  stabilized  nuclear  cloud. 
The  larger  mass  in  submicron  particles  in  the  nu¬ 
clear  cloud  can  result,  for  example,  from  fireball 
vapoiiration  and  recondensation  processes  or 
from  efficient  pulverization  of  soil  aggregates  by 
the  blast. 


Since  the  nuclear  cloud  is  likely  to  be  much  thick¬ 
er  and  the  PSD  different,  we  calculate  the  column 
density  M  as  follows; 


..d 


M(d)  =  -  U 


R(g  d^ 


(11) 


where 


<*2 


R  =  - 


sM  , 

dt 


top 


3 

3z 


't>(zl  dz 


(12) 


Here  U  is  the  horizontal  transport  speed,  iMz)  is 
the  vertical  mass  flux  crossing  altitude  z.  and  zjop 
is  the  maximum  plume  altitude.  Horizontal 
spreading  and  vertical  wind  shear  is  neglected. 
For  a  given  PSD,  i.e,.  dm/d(log  r),  the  vertical 
mass  flux  is: 


Thus  R  is  equivalent  to  the  mass  flux  at  the  sur¬ 
face. 


We  use  Eqs.  (11),  (12),  and  1 1 3)  to  compute  M 
as  a  function  of  d  for  several  hypothetical  nuclear 
PSDs.  Initial  dust  concentrations  and  PSDs  are 
specified  and  assumed  to  be  independent  of  alti¬ 
tude  throughout  the  troposphere  (z  <  10  kmi ;  no 
stratospheric  loading  is  assumed.  The  tropo¬ 
sphere  is  divided  into  ten  layers  i  Az  =  i  km i .  Par¬ 
ticles  in  seven  radius  bins  [r^in  =  0.01  um.  r^nax 
=  32  um.  Adog  r)  =  0.5]  fall  through  the  inter¬ 
faces  between  these  layers  at  a  sedimentation 
speed  V,  (proportional  to  r'^2)  given  by  Stokes  s 
law.  Separate  mass  balances  are  kept  for  each 
layer.  As  each  layer  “feeds"  dust  into  the  neigh¬ 
boring  layer  below  it  (and  receives  dust  from  the 
layer  above),  the  spectral  mass  concentration 
dm/ddog  r)  changes.  The  net  change  in  column 
density  for  each  layer  over  a  small  time  interval 
At  is  <D(z)  At.  Integrating  this  quantity  over  the 
entire  dust  column  leads  to  a  time-dependent  es¬ 
timate  of  R.  This  time  dependence  is  converted 
to  a  downstream  spatial  dependence  by  the 
transformation  Ad  =  UAt.  where  U  is  the  horizon¬ 
tal  transport  speed.  We  then  compute  Mfd)  is 
from  Eq.  (II).  We  assume  U  =  20m.'s,  consistent 
with  Saharan  dust  transport. 


Four  initial  PSDs  were  tested.  Three  are  based 
on  a  simple  power  law: 


ddog  r; 


max 

- 

<D(Z)  = 

V,(9.z)^ 

d9 

- 

(13) 


where  Vj(0,z)  is  the  sedimentation  velocity  (a 
function  of  particle  radius  r).  9  =  log(r).  and  ©max- 
9min  correspond  to  the  upper  and  lower  radii 
bounds.  Substituting  Eq.(12)  into  Eq.(11)  and 
using  the  condition  <t>(Ztop)  =  0  gives 


(•  max 


R  = 


9  , 
min 


V,(9) ^ 
d0 


de 


surface 


(14) 


where  A  and  0  are  constants  and  r  <  32  um. 

For  3  =1 .0,  dm/d(log  r)  resembles  a  typical  Sa¬ 
haran  source  region  PSD  (see,  for  example,  the 
upper  curve  in  Fig,  7),  The  submicron  fraction  for 
this  distribution  is  3. 1  percent  by  mass.  For  com¬ 
parison,  we  also  tested  two  cases  with  the  same 
initial  mass  concentration  (700  pg/mT),  but  a  fin¬ 
er  distribution  (3  =  0.5;  submicron  fraction  =16.3 
percent)  and  a  coarser  distribution  (3  =  1 .5;  sub¬ 
micron  fraction  =  0.6  percent).  These  PSDs  are 
plotted  in  Fig.  12.  A  least-square  fit  to  the  PSDs 
observed  in  the  Johnie  Boy  stabilized  cloud  [Na¬ 
thans,  Thews,  and  Russell,  1970]  was  also  used 
as  an  initial  condition.  For  r^g^  =  324m.  the  Joh¬ 
nie  Boy  submicron  fraction  is  1 9  percent  by  mass 
(normalized  such  that  the  initial  mass  concentra¬ 
tion  is  700  4g/m3) . 
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Figure  12.  Three  hypoihelical  power  law  PSDs  for  low  altitude  micleai  buists  and  least  squares 
subsurface  burst  PSD  (mass  concentration  is  700  mj  iiD) 


The  normalized  column  density  for  all  four  cases 
was  computed  for  a  total  transport  time  of  five 
days.  For  U  =  20  m/s.  this  is  equivalent  to  d  <8 
8600  km.  The  results  are  shown  in  Fig.  13.  In  all 
cases  the  exponential  dependence  of  M  on  d  is 
qualitatively  similar  to  the  Saharan  case  (Fig. 5). 
However,  it  is  dear  from  Fig.  13  that  the  Johnie 
Boy  distribution  had  the  greatest  concentration 
of  micron-size  particles.  This  is  evident  in  the  rap¬ 
id  initial  decrease  in  column  density.  The  subse¬ 
quent  decay  is  much  slower— almost  as  if  the 
initial  PSD  was  bimodal.  Indeed,  the  Johnie  Boy 
distribution  was  more  consistent  with  a  buried  or 
surface  burst  than  a  near-surface  air  burst  PSD 
'see.  for  example.  Pinnick,  Fernandez,  and 
Hinds.  1983). 

There  is  a  strong  dependence  of  the  depletion 
rate  on  the  initial  PSD.  In  particular,  the  greater 
the  initial  concentration  of  large  (r  >  1  nm)  par¬ 
ticles,  the  greater  the  initial  depletion  rate.  On  the 
other  hand,  since  these  particles  are  rapidly  de¬ 
pleted.  the  late-time  dust  cloud  consists  largely 
of  submicron-size  particles.  Hence,  for  large  d, 
M  is  proportional  to  the  initial  submicron  fraction: 
at  d  =  8000  km,  M  (p  =  0.5)  is  nearly  five  times 
larger  than  M(3  =  1.5).  The  change  in  M  over 
most  of  the  lifetime  of  a  nuclear  cloud  subject 
only  to  dry  deposition  is  determined  primarily  by 
the  initial  concentration  of  submicron-size  par¬ 
ticles.  Larger  particles  have  only  a  transient  ef¬ 
fect  upon  the  mass  deposition  rate. 


Late-time  depletion  is  accelerated  by  coagulation 
and  precipitation.  These  effects  were  nnt  in¬ 
cluded  in  the  four  calculated  cases.  Rain— paiuc- 
ularly  over  the  North  Atlantic— would  accelerate 
dust  removal  from  the  troposphere.  Similarly,  co¬ 
agulation  would  accelerate  the  sedimentation 
rate.  Finally,  Eqs.  (11).  il2).  and  (13)  account 
for  neither  lateral  spreading  <  diffusion  i  of  the  nu¬ 
clear  cloud  nor  wind  shear.  Nevertheless,  the  re¬ 
sults  demonstrate  the  strong  sensitivity  of  mass 
removal  rates  to  changes  in  the  initial  PSD. 

RADIATIVE  EFFECTS. 

Saharan  plumes  are  highly  reflective  at  solar 
wavelengths,  and  absorption  is  rather  weak. 
These  characteristics  are  dictated  by  ( i )  the 
composition  of  the  dust  particles  (primanlv  sili¬ 
cates).  and  (2)  the  size  of  the  particles.  Nuclear 
dust  clouds  that  have  a  similar  composition  and 
size  distribution  are  likely  to  have  similar  radiative 
properties.  However,  dust  lofted  from  U.S,  and 
Soviet  silo  fields  would  have  a  heavier  organic 
content  than  Saharan  dust,  resulting  m  greater 
absorption  and  increased  cloud  optical  depths. 
As  a  consequence,  there  may  be  some  warming 
in  the  upper  troposphere  and  stratosphere.  This, 
in  turn,  would  tend  to  stabilize  the  atmosphere 
and  possibly  suppress  mid-latitude  storm  activity. 
A  similar  mechanism  has  been  observed  on  Mars 
during  major  dust  storms  [Haberie.  1986|. 
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e  13.  Computed  normalized  column  densities  as  function  of  dista  ice  downsiieam  (or  toui 
initial  PSDs. 


SECTION  5 


CONCLUSIONS 


We  find  that  the  mass  of  Saharan  dust  lofted  is 
comparable  to  about  one-third  the  injection  ex¬ 
pected  from  a  full-scale  counterforce  attack 
against  a  U.S.  Minuteman  missile  field.  The 
amount  of  dust  raised  is  highly  dependent  upon 
soil  condition  and  type  and  can  vary  significantly 
with  season.  Seasonal  changes  in  vegetation,  soil 
moisture,  and  the  presence  of  surface  crust  re¬ 
duce  the  scouring  efficiency  in  Saharan  source 
regions.  Similar  dependencies  are  expected  in 
U.S.  and  Soviet  silo  fields  where  cultivation, 
vegetation  growth  cycles,  ground  freeze,  and 
snow  cover  are  factors.  In  the  nuclear  case,  soils 
are  further  modified  by  thermal  and  blast  effects. 

Lofting  of  Saharan  dust  is  qualitatively  similar  to 
nuclear  sweepup.  in  both  cases,  strong  surface 
winds  are  required  to  scour  dust  from  the  surface 
and  a  buoyant  momentum  source  is  needed  to 
carry  the  dust  to  altitude.  The  strength  of  these 
mechanisms  is  several  orders  of  magnitude  larg¬ 
er  in  the  nuclear  case.  Consequently,  a  wider 
spectrum  of  particle  sizes  is  initially  lofted  and 
carried  to  greater  altitudes.  The  largest  particles 
are  rapidly  depleted:  therefore,  after  several 
hours,  particle-size  distributions  and  mass  con¬ 
centrations  in  nuclear  clouds  are  similar  to  those 
observed  in  Saharan  plumes. 

Saharan  dust  plumes  extend  4  to  6  km  in  altitude. 
Nuclear  clouds  may  reach  the  stratosphere.  De¬ 
spite  the  relatively  shallow  injection.  Saharan  dust 
is  transported  several  thousand  kilometers  down- 
stteam.  The  long-range  transport  is  aided  by  (1) 
the  presence  of  a  low-altitude  stable  layer  which 
isolates  the  dust  from  turbulent  surface  pro¬ 
cesses.  and  (2)  the  absence  of  wet  removal 
(precipitation  scavenging).  Mass  depletion  rates 
are  initially  "*200  kg/km2/h:  this  value  decays 
rapidly  downstream  as  the  most  massive  par¬ 
ticles  are  removed  by  sedimentation.  After  2000 


km,  roughly  70  percent  of  the  initial  dust  mass 
is  depleted.  At  greater  distances,  depletion  rates 
approach  2  kg/km^'h:  the  exact  value  is  propor¬ 
tional  to  the  initial  concentration  of  submicron-ra- 
dius  particles  and  boundary  layer  turbulence.  Nu¬ 
clear  clouds,  which  may  contain  a  large  mass  of 
very  fine  recondensed  particles,  would  retain 
more  of  their  initial  mass  at  late  times,  in  mioole 
latitudes,  where  nuclear  clouds  originate,  precipi¬ 
tation  scavenging  is  likely  to  dominate  the  mass 
depletion.  This  scavenging  is  enhanced  when  the 
cloud  contains  large  numbers  of  submicron-size 
clay  ice  nuclei  particles.  Saharan  plumes  thus 
represent  an  upper  limit  to  the  longevity  of  dust 
in  the  lower  troposphere. 

Far  from  the  source,  very  large  particles  ir  >  10 
nm)  are  sometimes  found  in  the  Saharan  plume. 
This  may  be  evidence  of  cloud-aided  coagula¬ 
tion.  Nuclear  clouds,  which  encounter  moisture 
clouds  more  often,  may  contain  a  greater  con¬ 
centration  of  such  ultragiant  particles,  presenting 
a  hazard  to  aircraft  and  other  flight  vehicles. 

Saharan  dust  has  a  very  low  organic  content, 
while  dust  from  likely  target  areas  would  have  a 
greater  organic  content.  Nuclear  c'ouds  are 
therefore  more  absorptive  (at  least  at  solar  wave¬ 
lengths)  than  Saharan  plumes.  W'e  also  find  tfnt 
nuclear  clouds  are  probably  less  reflective  than 
Saharan  plumes  due  to  greater  concentrations  of 
submicron-size  particles.  Optical  depths  in  the 
more  massive  Saharan  outbreaks  may  approxi¬ 
mate  those  expected  in  nuclear  clouds. 

Saharan  plumes  provide  an  ideal  natural  ana¬ 
logue  for  the  late-time  removal  of  dust  from  the 
lower  troposphere  in  the  absence  of  precipita¬ 
tion.  Removal  rates  as  a  function  of  downstream 
distance  are  well  determined.  Saharan  plumes 
are  therefore  appropriate  test  cases  lor  nuclear 
dust  transport  models. 
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